Majorana fermions in spin-singlet nodal superconductors with coexisting non-collinear 

magnetic order 
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Different realizations of Majorana fermions in solid state materials have attracted a great deal 
of interests recently in connection to topological order and quantum information processing. We 
propose a novel way to create Majorana fermions in superconductors. We show that an incipient 
non-collinear magnetic order turns a spin-singlet superconductor with nodes into a topological su- 
perconductor with a stable Majorana bound state in the vortex core or on the edge. We argue 
that such an exotic non-Abelian phase can be realized in extended t-J models on different lattices 
with two explicit examples: one on the triangular lattice and the other on the square lattice. Our 
proposal provides a new avenue for the search of Majorana fermions in correlated electron materials 
where nodal superconductivity and magnetism are two common caricatures. 

PACS numbers: 71.27. +a, 74.25.Ha, 74.81. Bd 



A Majorana fermion is an electrically neutral fermion 
whose antiparticle is itself. It was proposed to exist in 
high-energy physics such as neutrino physics and dark 
matter [1]. In recent years, Majorana fermions have at- 
tracted growing attention in condensed matter physics 
[2-4]. Specifically, Majorana fermions can be realized as 
zero-energy bound states in the vortex core or on the 
edge of certain two-dimensional superconductors. In- 
stead of the usual Bose or Fermi statistics, these vor- 
tices obey non-Abelian statistics [5-10] as a manifes- 
tation of topological order [11, 12]. Due to this re- 
markable feature, Majorana bound states (MBSs) can 
be utilized for topologically-protected qubits in fault- 
tolerant quantum computation [2, 13, 14]. Several sys- 
tems have been proposed to realize MBSs, such as 
even-denominator fractional quantum Hall states [5, 6, 
8, 15], p + ip superconductors [8- 10] and supcrfiuids 
[16, 17], superconductor-strong-topological-insulator in- 
terfaces [18-20], s-wavc Rashba superconductor [21-23] 
and spin-orbit-coupled nodal superconductors [24]. 

In this work, we present a novel realization of MBSs in 
spin-singlet superconductors with nodal excitations. We 
show that when a coexisting non-collinear magnetic or- 
der with a wavevector connecting two nodes of opposite 
momenta develops, there will be one MBS in each vor- 
tex core and on the edge of such topological supercon- 
ductors. We demonstrate our proposal with two explicit 
examples. The first one is a nodal d+id superconductor 
[25] coexisting with 1x3 (or 3x3) coplanar magnetic 
order on the triangular lattice. We argue that this state 
is likely to be realized in a doped t- Ji model on the tri- 
angular lattice and is relevant for a material: supercon- 
ducting sodium cobaltate Na^CoO^ • 2/H2O near doping 
concentration x = 1/3 [25, 26]. The second example is 
a d x 2_ y 2 superconductor with coexisting Q = (Qo.Qq) 
coplanar magnetic order on the square lattice which may 
be realized in a doped t- J\- J2- J3 model on the square lat- 
tice. Experimentally the phase diagram of a wide range 



of strongly correlated materials, from high-T c cuprates 
to heavy- fermion compounds, exhibits d-wave supercon- 
ductivity [27, 28] in proximity to [29] or coexisting with 
[28, 30] magnetic orders. These results suggest that Ma- 
jorana fermions may exist in correlated electron materials 
with magnetic frustration and nodal superconductivity. 

We begin with a general discussion of our theory. The 
low-energy excitations of a nodal superconductor are 
massless Dirac fermions with linear dispersion around 
the nodes. Our basic idea is to create a topological 
superconductor with a single MBS in the vortex core 
by adding a proper mass term to the nodal fermions. 
Let us consider a spin-singlet superconductor with n 
pairs of isolated nodes located at crystal momenta ±q^, 
£ = 1, . . . n. Expanding around the nodes, the low-energy 
BCS Hamiltonian describing the quasi-particle excita- 
tions has a generic form in the Nambu basis "I^k = 
(cq.+k/bclq^k^c-q.+k^-c^^) 11 for each pair of 
nodes at opposite momenta: 



"H off = J^^k^k^k, H ek = (n £k -t)ct z , (1) 



(ReA qf+k , -ImA qf+k , £ q(! +k) 
k a vi + 0(\kf), 



where Ak is the pairing gap function and £k is the ki- 
netic energy f and a are Pauli matrices operating in the 
particle-hole (Nambu) and spin sectors respectively. A 
particle- hole unitary rotation U = exp[i0f • f\ turns 
into U* HnJJ = {k\T x + &2' r y) (J z where /c 12 arc linearly- 
independent combinations of k x and k y [51]. 

Let's focus on a single ^-th pair of nodes at ±qj. 
Clearly, a gap will be generated by adding a generic mass 
term To(a x KeAl — a y lmM), where M is a complex order 
parameter, in the Hamiltonian for this pair of nodes, 

tfH'nJJ = (kiT x + k 2 Ty)a z +T (a x ReM -a y lmM). (2) 
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This is nothing but the effective Hamiltonian for s-wave 
superconductivity on the surface of a 3D topological in- 
sulator with M playing the role of the superconducting 
order parameter [18]. It is known that the latter contains 
a single MBS in its vortex core. In the present context of 
singlet nodal superconductors, the physical origin of the 
local order M turns out to be a non-collinear (co-planar) 
magnetic order described by 

H cp = [M(S* + iS»)e a *>- r + h.c] 

r 

= X! MC L+k+qo,t C ^+ k - c l"-i + /l ' C -' ( 3 ) 

where = ^2 a g cl a a^ ^c rjQ ,/2, a = x,y,z are the elec- 
tron spin operators at lattice site r. When the ordering 
wave vector 2qo connects the pair of nodes located at 
±q^, the magnetic scattering generates precisely the mass 
term in Eq. (2). Remarkably, this non-collinear mag- 
netic order drives the nodal superconductor into a non- 
Abelian topological state since, among the two (even and 
odd) combinations c kj£ ( ) = ^(c qf+k ,t ± c -q £ +k,j.), the 
odd-combination is driven into the weak-pairing phase 
[8] by the mass gap, while the even-combination enters 
a strong-pairing phase. The situation is analogous to a 
doubled-layer v — 1/2 fractional quantum Hall system, 
where an Abelian (331) state can be driven to a non- 
Abelian pfaffian state by a tunneling term between the 
layers [8, 31, 32]. In contrast to the pair of nodes with 
opposite momenta connected by the magnetic ordering 
vector 2qo, the other nodal fermions are spin- flip scat- 
tered either to finite energy states away from the Fermi 
level or, in certain special cases, to different nodes not 
connected by the pairing part. As a result, they either 
remain gapless or be gapped out in the magnetic sector. 
We note that the existence of a single MBS will not be 
affected by the presence of other nodal excitations [24]. 

On general grounds, the non-collinear magnetic order 
can be realized in a frustrated system from the residual 
spin-spin interactions between the nodal fermions and 
coexists with superconductivity. We stress that it is cru- 
cial to require the magnetic order to be non-collinear: a 
collinear spin order such as 

% c i = 2m^5' r a; cos(2q -r) 

r 

= S5Z mC q*+k+q 0l <rC«U+k-qo 1 ff + '*-C. ) (4) 
£k a 

will not only drive the Nambu pair (c qo +k,t, c l qo -k {)■> 

but also (c_ qo +k.t: c qo - k 4.) m t° the weak-pairing phase. 
This amounts to two copies of weak-pairing p + ip super- 
conductors. Thus, there will be two MBSs in the vortex 
core and two counter-propagating Majorana modes on 
the edge. The two branches will scatter with each other 
and open up a gap in the edge spectrum. As a result, 




FIG. 1: (color online) (a) The 1st BZ of the triangular lattice 
and the 6 nodes (Ni) of the 2nd NN d + id pairing gap func- 
tion. The normal state hole FS (dashed circle) crosses the 
gap nodes at doping x c for the nodal chiral superconductor. 
The arrow indicates the momentum of the 1x3 magnetic or- 
der with the magnetic zone shown by dashed parallelogram, 
(b) The spin configuration of the 1x3 non-collinear magnetic 
order. Si^ are two primitive lattice vectors. The reciprocal 
vectors are 61,2 with the convention 3i ■ bj = 5i t j. 



there will be no MBS in a nodal superconductor with 
collinear magnetic order. 

In the presence of the non-collinear magnetic order (3), 
both inversion symmetry and spin rotational symmetry 
are broken and thus the spin-singlet pairing can mix with 
spin triplet pairing. When the triplet pairing amplitude 
is small compared with |M|, the system will stay in the 
gapped non-Abclian topological phase. In the opposite 
limit, a dominant one-component chiral triplet pairing 
state is well known to be in the non-abelian weak-pairing 
phase [8, 9]. Therefore we expect that the non- Abelian 
topological superconductor to be stable against the mix- 
ing between singlet and triplet pairing. In the following 
we shall demonstrate the above general discussions by 
two specific examples on the triangular and the square 
lattices. 

We start with the nodal chiral superconductor with 
coplanar 1x3 magnetic order on the triangular lattice. 
Nodal chiral superconductors were proposed in the con- 
text of the sodium cobaltate superconductivity [25] . The 
sodium cobaltate (Na^OoC^) is a layered, strongly cor- 
related triangular lattice electron system which becomes 
a superconductor upon water-intercalation at low tem- 
peratures [33]. Recent NMR measurements finds strong 
evidence for singlet pairing [26, 34, 35] with nodal excita- 
tions at a critical doping x c ~ 0.26 [35], which can be ex- 
plained by the nodal d+ id pairing state [25]. Specifically, 
2nd nearest neighbor (NN) d+ id pairing was argued to 
be the dominant pairing channel on the electron doped 
triangular lattice where the complex gap function has 
6 isolated zeros inside the 1st Brillouin zone (BZ). The 
Fermi surface (FS) crosses these nodes at a critical dop- 
ing concentration x c producing 6 Dirac points as shown 
in Fig. 1(a). The superconducting states at x < x c and 
x > x c arc separable by a topological phase transition. 
We therefore consider a simple effective pairing model 
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described by the Hamiltonian, 

U = X] £k4 CT c ka + ^(A k Ckt-C-14 - h.c), (5) 
kcr k 

where £k is the dispersion for electrons hopping with 
amplitudes {ti,t 2 ,t^) = (—202,35,29) meV between the 
first 3 nearest-neighbors (NNs) [25]. A k = 2A 2 [cos(A: 1 - 
k 2 ) + e i27r / 3 cos(2A:i + k 2 ) + e i4 * /3 cos(fci + 2k 2 )] is the 
2nd NN d + id pairing gap function in the basis k = 
kibi + k 2 b 2 shown in Fig. 1(a) where the 6 Dirac nodes 
(Ni, i = 1, • • • ,6) are located at {ki,k 2 ) = ±(2tt/3, 0), 
±(0,2tt/3), and ±(2tt/3, -2tt/3). Such a d + id su- 
perconductor exhibits quantized spin Hall conductance 
[8] associated with the winding number W of the unit 

(2) (2) 

vector fik = (RcA k , — ImA k , £k)/-E k where £k = 

+ | A k 2 ^ | 2 . When the FS lies inside the Dirac points 
(x > x c ), W = — 2 and there are two counter-clockwise- 
propagating chiral fermions on the edge; each is charge 
neutral but carries spin h/2 [25]. When the FS encloses 
the six gap nodes (x < x c ), W = 4 and there are four 
spin-carrying chiral fermions on the edge. 

Remarkably, a non-collinear magnetic order described 
by Hep in Eq. (3), with the vector 2qo pointing from 
iVj+3 to Ni (Fig. la), produces a mass gap for the nodal 
fermions as discussed above. The magnetic order corre- 
sponds to the 1x3 coplanar pattern shown in Fig. 1(b). 
The nodal chiral superconductor is thus turned into a 
non-Abelian topological superconductor with a winding 
number W = +1, which is similar to a spinless p+ip su- 
perconductor [8] and supports a single MBS in the vortex 
core and on the sample edge. To demonstrate the latter, 
we calculate explicitly the edge spectrum of such a non- 
collinear magnetic superconductor described by % + % cp 
with two parallel edges along the (1,1) direction. The re- 
sults arc shown in Fig. 2 for A 2 = 150 meV and magnetic 
mass M = 100 meV. The bulk excitations are completely 
gapped, due to the fact that the magnetic scattering wave 
vector 2qo not only connects the Nambu pair (N^,N 2 ) 
but also (N3,Ni) and (Ni,N e ) in the magnetic sector. 
The in-gap edge states are separated from the bulk con- 
tinuum. The central result is that there is a single branch 
of gapless edge mode: a single MBS crossing k = local- 
ized at the two edges respectively. We verified that, as 
discussed in the generic case above, a collinear magnetic 
order described by H c i in Eq. (4) will not create a MBS 
in the edge spectrum, although the bulk excitations are 
fully gapped (see supplemental materials). 

A natural question rises: how does the 1x3 non- 
collinear magnetic order arise microscopically on the tri- 
angular lattice? To this end, we consider the 2nd NN 
antiferromagnetic Heisenberg model, i.e. the J 2 model, 
on the triangular lattice. The classical ground states of 
the J 2 model is well-known to have 3x3 non-collinear 
magnetic order [36, 37]: on each of the three sublattices 
connected by 2nd-NN bonds the spins exhibit 120 degree 




k: momentum along b-2 direction 

FIG. 2: (color online) Edge state spectra as a function of mo- 
mentum k along &2 direction in the magnetic zone as shown 
in Fig. 1. The in-gap edge states are separated from the bulk 
continuum. Blue and red lines correspond to the edge states 
localized on the two separate edges respectively. Note that 
there is a single counter-clockwise-propagating zero-energy 
MBS at k = 0. 



coplanar order. There is a large ground state degeneracy 
at the classical level from the relative spin orientations on 
the three sublattices. Such a 3 x 3 noncollinear magnetic 
order as a classical ground state of J 2 model already in- 
duces a single MBS crossing k = in the edge spectrum, 
as shown in the supplemental materials. Quantum fluc- 
tuations would further lift this degeneracy at classical 
level through the order-due-to-disorder mechanism [38- 
40], such that the true quantum ground state has 1x3 
order. We have carried out a Schwinger-boson large-.!? ex- 
pansion study (see supplementary materials) of the spin- 
s' Heisenberg J 2 model and found that when S is larger 
than a critical value S c « 0.17 (such as for spin-1/2), the 
system develops the 1x3 non-collinear order as shown 
in Fig. 1(b). This suggests that if the residual interac- 
tion between the nodal fermions is dominated by the 2nd 
NN Heisenberg exchange J 2 , then the 1x3 non-collinear 
magnetic order is likely to develop with the wave vector 
2qo connecting the gap nodes with opposite momenta as 
shown in Fig. 1. Even more intriguingly, to the extent 
that J 2 would favor a 2nd NN d + id resonance valence 
bond pairing state, it is likely that both the nodal chiral 
superconductor and the non-collinear magnetic order can 
emerge from the same exchange interaction in a doped t- 
J 2 model. 

We turn to the 2nd example of a more familiar nodal 
d x 2_ y 2 superconductor with (Qo,Qo) coplanar magnetic 
order on the square lattice. The Hamiltonian has the 
same form as T~i in Eq. (5), but with the dispersion 
£k = — 2i[cos(fci + k 2 ) + cos(fci — k 2 )] — fi for NN hop- 
ping t, and the NN d x 2_ y 2 pairing gap function Ak = 
2Ai [cos(fci + k 2 ) — cos(fci — k 2 )~\ . The momentum is de- 
fined as k = fci&i + k 2 k 2 as shown in Fig. 3, where four 
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(b) 



FIG. 3: (color online) (a) The 1st BZ of the square lattice. 
Black circles denote the four nodal points (Ni) where the 
FS intersects the nodal lines of the NN d x 2_ y 2 pairing gap 
function. The non-collinear magnetic order with momentum 
(7r/2, 7r/2) (shown by the arrow) turns the node-structure into 
the 6 red diamonds inside the magnetic zone (red rectangle), 
(b) The spin configuration of the (n/2,n/2) non-collinear 
magnetic order. ai,2 represent two primitive vectors. The 
reciprocal vectors are 61,2 with Si ■ bj = 5ij. Along E2 direc- 
tion the unit cell is doubled due to magnetic order. 



nodal points arise from the crossing of the FS with the 
nodal lines of the gap function: N\$ : (ki — 0, k% — ±<Zo) 
and iV2,4 : (k% = ±go,fc2 = 0) with go = arccos(— ^). A 
non-collinear magnetic order described by H cp in Eq. (3) 
with ordering momentum Qo = 2qo = 2go&2 gaps out 
the (cjv S) -f-, c' Ni 1 ) branch and creates a single MBS. Fig. 3 

shows a specific example with /i = — 2\/2t and go = 7r/4, 
together with the spin configurations. In this case, the 
commensurate magnetic order cannot gap out all nodal 
excitations since the Hamiltonian is still invariant un- 
der time reversal followed by a lattice translation [41]. 
As shown in Fig. 3(a), the non-collinear magnetic or- 
der splits the original four spin-degenerate nodal points 
(black circles) into 6 non-degenerate nodes (red dia- 
monds). The vanished pair of nodes is gapped out by 
the magnetic mass (3) and enters the weak-pairing phase. 
The calculated edge spectra along the (1,1) direction is 
shown in Fig. 4 near the magnetic zone boundary. It 
confirms the existence of zero energy MBS localized on 
the parallel edges. Similar results are obtained at other 
commensurate values of go such as go = 7r /3. Since the 
gapless bulk excitations are located at different momenta, 
the MBS near ki = ir is expected to be stable against im- 
purities and the mixing with bulk excitations [24]. For a 
generic doping, go is incommensurate with the lattice and 
a corresponding incommensurate non-collinear magnetic 
order could open up a full gap for the bulk excitations 
and leave the edge states as the only low-energy modes. 

A surprising and remarkable feature seen in Fig. 4 is 
that the Majorana mode on the edge is dispersionless, 
i.e. it is localized and does not possess a chirality. It 
turns out that the chirality of the Majorana mode on the 
edge crucially depends on the symmetry of the pairing 
order parameter Ak in (1). The chirality of the edge 
mode associated with the ^-th pair of gap nodes can be 
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k2'. momentum along (1,1) direction 

FIG. 4: (color online) Edge spectra as a function of momen- 
tum along the (1, 1) direction of the magnetic zone for 
a NN d x 2_ y 2 superconductor with (tt/2,tt/2) non-collinear 
magnetic order. The spectrum is zoomed in around ki — it 
where a flat band of zero energy states are localized on the 
two edges (blue and magenta), corresponding to the MBS. 
Red dispersing lines are the bulk states and the two nodal 
points in the bulk excitations are located at k^ ~ 7r ± 0.04. 



where v x y are the 



defined as C g = sign[z • (v x x flj)] 
velocities defined in (1). It is straightforward to verify 
that this definition is invariant under a global U(l) rota- 
tion Ck,cr — > e ie Ck,cr. For a singlet nodal superconductor 
that preserves the time-reversal symmetry, which is the 
case for the d x 2_ y 2 superconductor, the chirality of the 
edge mode C = and hence the non-dispersing Majo- 
rana edge mode in FIG. 4. However in the presence of 
the magnetic order (3), a time reversal symmetry break- 
ing pairing term can be induced. We verified that adding 
a small imaginary part in Ai indeed generates a single 
MBS on the edge that disperses across k2 = tt with cer- 
tain chirality (see supplemental material). 

It is possible to realize this non-collinear magnetic or- 
der in the microscopic Heisenberg J1-J2-J3 model on the 
square lattice. The classical ground state of the J1-J2-J3 
model is a non-collincar one with momenta Q = Qo^i,2 
where cos(Q ) = - J 1 /(2J 2 + 4J 3 ) for 4J 3 + 2J 2 > Ji and 
•h > ^2/2 [42, 43]. There are numerical evidence that the 
latter survives in the quantum 5 = 1/2 Heisenberg Ji- 
J2-J3 model in a wide parameter range [44, 45]. There- 
fore we expect that the non-Abelian phase of a nodal 
d-w&ve superconductor coexisting with Qo = Qo&i,2 non- 
collinear magnetic order may be realized in certain pa- 
rameter regime of the doped t- J1-J2-J3 model. 

In summary, we proposed a new type of non-Abelian 
topological superconductors. They emerge when spin- 
singlet superconductors with isolated nodes coexist with 
non-collinear magnetic order with the ordering vector 
connecting the nodes at opposite momenta, due to the 
exchange interaction between the nodal fermions. Majo- 
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rana fermions arise in the vortex core and on the edge 
of such magnetic superconductors. Since magnetism and 
unconventional superconductivity are common features 
of strong correlation, our findings suggest searching for 
the non-Abelian MBSs in correlated electronic materials 
with magnetic frustration and nodal superconductivity. 

We thank Sen Zhou for useful discussions. This work 
is supported in part by NSF DMR-0704545 and DOE 
DE-FG02-99ER45747. 
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A. Noncollinear magnetic order of J2 Heisenberg 
model on triangular lattice 

Due to geometric frustration, the triangular lattice 
antifcrromagnet has long been considered as a promis- 
ing candidate for realizing the spin liquid state [46]. 
These spin disordered states are different from conven- 
tional symmetry-breaking phases and its low-energy exci- 
tations are deconfincd fractionalized quasiparticles (such 
as bosonic/fermionic spinons which carry spin but not 
electric charge). Strong quantum fluctuations at low tem- 
peratures may stabilize such disordered states; whereas 
in the large- S classical limit, various symmetry breaking 
magnetic ordered states develop with their low-energy 
physics dominated by Goldstone modes associated with 
the broken symmetry, i.e. spin waves or magnons. A 
well-known example is the zero-flux state introduced by 
Sachdev [47] in the Schwinger-boson representation. As 
the spin S (or the occupation number of bosonic spinon 
per site) increases to a critical value [48], the spinons at 
the zone corner condense, leading to the 120 degree or- 
dered state as predicted for the classical Heisenberg Ji 
model on the triangular lattice. In general a continuous 
quantum phase transition between a disordered spin liq- 
uid and a non-collinear magnetically ordered phase can 
be described by Bose condensation of spinons [49]. 

Our strategy here is to use the Schwinger-boson repre- 
sentation to study the Z2 spin liquids for the 2nd nearest 
neighbor (NN) J 2 Heisenberg model on the triangular 
lattice. We then track down the pattern of the non- 
collinear magnetic order by studying the Bose conden- 
sation of spinons in the neighborhood of the spin liquid 
state using the Schwinger boson mean field theory. 



Schwinger-boson representation of Z2 spin liquids in 
triangular lattice J2 model 

In the Schwinger-boson representation the spin opera- 
tors S r on lattice site r are written in terms of Schwinger 



bosons b ra 

S r = \ b t,*v<x,0 b r,P ( 6 ) 

where a denotes the three Pauli matrices. The spin value 
S is fixed by the constraint 

b UK« = « = 25, V r. (7) 

Here k = 2S is a parameter measuring the strength of 
quantum fluctuations in the system and k — > +00 corre- 
sponds to the classical limit. Consider two sites i and j, 
there are only two bond variables that preserve the spin 
rotational symmetry 

Bij =lEa b i,oPi,<* 

where e Q( g is the rank-2 fully anti-symmetric tensor. 
Thus, a generic Heisenberg Hamiltonian can be written 
as 

H = J2 JijSi -S, X M-MAv + 440- (8) 

At the mean-field level, the pairing amplitudes Aij = 
(Aij) = —Aji and the hopping amplitudes By = (Bij) = 
Bji where Aij,B%j are complex variational parameters. 
Enforcing the constraint (7) globally by a chemical po- 
tential fx, we obtain the mean-field Hamiltonian 

h mf = Eij Jiji-AijAij + B i:j B^ + h.c.) (9) 

+ Ei,j Jiid^l 2 - \Bij\ 2 ) -mEt(»U« - K ) 

Different types of the spatial "patterns" of the am- 
plitudes {A^, Bij} correspond to different universality 
classes of the Z2 spin liquids. They are characterized 
by different symmetry-protected topological orders and 
classified by the projective symmetry group (PSG) [50]. 
In Ref. 48, all different Z 2 spin liquids preserving the lat- 
tice symmetry on the triangular lattice are classified in 
the Schwinger-boson representation. There are in total 
8 different PSGs on an isotropic triangular lattice, corre- 
sponding to 8 different universality classes of the Z2 spin 
liquids. 

Let's focus on those spin liquid states that are possible 
ground states in the 2nd NN J 2 Heisenberg model on 
an isotropic triangular lattice. We require the 2nd NN 
mean-field pairing amplitudes to be nonzero: Aij 7^ 
for ((ij)). This constraint excludes 6 of the 8 possible 
Z2 PSGs. Therefore we have only two different Z2 spin 
liquids with nonzero 2nd NN Ay 's: they are labeled as 
the 7r-flux state (following the name given in Ref. 48) and 
0-flux-2 state (in order to distinguish from the zero-flux 
state in Ref. 48). The 0-flux-2 state corresponds to the 



FIG. 5: (color online) The sign structure of 2nd NN pairing 
amplitudes Ay = h(bi-fbj± ~~ f° r the two different Z2 

spin liquids: (a) 0-flux-2 state and (b) 7r-flux state. Specifi- 
cally, writing Aij = ±iA2 with A2 a constant, then the plus 
sign is assigned if i — > j is along the direction of the triple 
arrows. ai,2 are two primitive lattice vectors. Note that the 
7r-flux state breaks translation symmetry along the 0,2 direc- 
tion and doubles the unit cell. 



*= 2S 



Z2 PSG solution p\ = 0, P2 = P3 = 1 in the notation of 
Ref. 48. Among them the 7r-flux state doesn't allow 2nd 
NN hopping terms, i.e. Bij = for all The 0-flux-2 

state, however, does allow uniform real 2nd NN hopping, 
i.e. Bij eB 2 £1 for all ((ij))- For both phases, the 2nd 
NN pairing amplitudes can be made purely imaginary by 
a gauge choice, with their spatial sign structures shown 
in Fig. 5. 

We carry out self-consistent calculations of the mean- 
field Hamiltonian (9) for both the 7r-flux and the 0-flux-2 
states. For simplicity we only include the 2nd NN pairing 
amplitudes for the O-flux-2 state. Thus, strictly speaking 
the obtained ground state energy for the O-flux-2 state 
gives an upper bound of the actual energy. The com- 
parison of the ground state energy (or zero-temperature 
free energy) for the two states are shown in Fig. 6, which 
shows that the O-flux-2 state has a lower energy than the 
7r-fiux state. A nonzero 2nd NN hopping amplitudes B2 
would lower the energy of the O-flux-2 state even further. 
Thus the lowest energy Z2 spin liquid state in the trian- 
gular lattice J2 model should be the O-flux-2 state which, 
as we will show below, is stable against spin order when 
the quantum parameter is smaller than n c w 0.34 (see 
Fig. 6). 



FIG. 6: (color online) The mean-field ground state energy per 
site E ga /{N S J2) = — 3| | 2 of the two Z2 spin liquid states in 
the J2 Heisenberg model on the triangular lattice. Blue aster- 
isks denote O-flux-2 state and red circles denote 7r-flux state. 
Apparently O-flux-2 state has lower energy than its competi- 
tor 7r-flux state and should be the ground state of J2 model 
for k < 0.34. When the quantum parameter n — 2S increases 
and becomes larger than critical values k c (O-flux-2) ~ 0.34 
and fv c (7r-flux) « 0.75, the spin system goes through a phase 
transition from spin liquid phases to magnetic ordered phases. 




FIG. 7: (color online) The spin configuration of the 1x3 
coplanar magnetic order in a S — 1/2 Heisenberg J2 model 
on the triangular lattice, obtained from spinon condensation 
in the O-flux-2 state. 



ten as 



Spinon condensation and the magnetic order in 
triangular lattice J2 Heisenberg model 



Now we show how to extract the magnetic ordering 
pattern in the O-flux-2 state when k becomes large. We 
choose bx,2 to be the reciprocal lattice vector which sat- 
isfies di ■ bj = Sij, i,j = 1, 2 where ai^ are the primitive 
vectors shown in Fig. 5. Expressing the momentum as 
k = fci&i + ^262, the mean-field Hamiltonian can be writ- 



H 



MF 



k4. 



+N S //(k + 1)+3J 2 |A 2 |' 
where N s is the number of lattice sites and 

(2) 



Ok.f 

6 -k4 



■9k 



(10) 



■ sin(fci — fc 2 ) — sin(/ci + 2^2) + sin(2fci + fc 2 ) 



7 



is the 2nd NN pairing structure factor corresponding to 
Fig. 5(a). The spinon dispersion, obtained from Eq. (10), 



Edge spectrum: d+id-wave SC coexisting with 3x3 SDW 



(2) I 



(11) 



has 6 minima at (fci,^) = ±(27r/3,0), ±(0, 27r/3) and 
±(27r/3, — 27r/3) which coincide with the 6 gap nodes 
{Ni} of the 2nd NN d + id pairing gap function [52]. 
We denote the momenta at points N{ as k;. As the 
quantum parameter k increases, the chemical poten- 
tial /i decreases and the minima of the spinon bands 
min{Ak} = \] l± 2 — 27| J2A2P/4 approach zero from the 
positive side. When fi = 3\/3>./2^l2/2 the band minima 
touch zero, the spinons at the 6 momentum points con- 
dense into 



1,2,3,4,5,6, (12) 



where Cj are complex numbers. The spatial configuration 
of the spinons is then given by 



" <6k 4 ,t> " 




1 




= Ci 


_(-!)*_ 



—x 



oj 



Etre 1 



oj 



—y 



03^3 



x-y 



(13) 



where r = xd\ + ya^ is the position vector of a site and 
lo = exp(i27r/3). In the above equation, 01,2,3 are three 
positive constants and £^1,2,3 are three SU(2) matrices 
with 



iUi 



c 

—c* 



j+3 



Ci+3 
C* 



1,2,3. 



The condensed spinons give rise to magnetic order and 
the corresponding spin configuration is 



(14) 



1 u 

,-2x 



Tr 



Tr- 



1 



U 



x Wi} 



uj- 2 y 1 

1 oj 2 (v-x) 



u. 



pu 3 }. 



In general, when all three ^'s are nonzero, the mag- 
netic order has the 3x3 spatial pattern on the trian- 
gular lattice. More precisely, there are three sublattices 
connected by the 2nd-NN J2 bonds: on every sublatticc 
the spins form coplanar 120 degree order as in the NN 
J\ Heisenberg model. When only one <f>i is nonzero, the 
corresponding spin configuration is the stripe-like 1x3 
coplanar ordered shown in Fig. 7. 

The classical ground state of the J1-J2 Heisenberg 
model with a = J2/J1 > 1 is known [36, 37] to have non- 
collinear magnetic order on the triangular lattice with 




k: momentum along ^2 direction 

FIG. 8: (color online) Edge spectrum of 2nd-NN d + id su- 
perconductor coexisting with an arbitrary 3x3 non-collinear 
magnetic order (14). Notations are the same as in Fig. 2 of the 
main text. Among the in-gap edge states separated from the 
bulk continuum, blue lines represent edge states on one edge, 
while red lines represent edge states on the other edge. The 
edge is along (1,1) direction. There is one single clockwise- 
propagating MBS crossing k = in the magnetic zone. The 
data are obtained for A2 = 150 meV and magnetic mass 
M — 100 meV for all three ordering vectors. 



the following momenta 
k = (fci,fc 2 



±(Q ,0), ±(0,Q ) and ±(Q ,-Qo), 
Qo = arccos(-i - i). (15) 



In the limit a — > +00, they exactly correspond to the 
magnetic order (14) obtained by condensing spinons in 
the 0-flux-2 spin liquid state. These classical solutions 
have a large ground state degeneracy, due to the relative 
spin orientations between different 2nd-NN sublattices as 
mentioned earlier. When quantum fluctuations are con- 
sidered, the degeneracy at the classical level will be lifted 
through the "order due to disorder" mechanism [38, 39]. 
Therefore the ground state of the quantum Heisenberg 
J2 model is expected [40] to exhibit magnetic order with 
one ordering wavevector among those in (15). In other 
words, only one of the {ci, i = 1, • • • , 6} in (13) is nonzero 
and the corresponding spin configuration (14) is the 1x3 
noncollinear order shown in Fig. 7. 



B. Majorana fermions with 3x3 noncollinear 
magnetic order 

It turns out that even for an arbitrary classical ground 
state (14) with a 3x 3 non-collinear magnetic order, when 
it coexists with the 2nd NN d + id superconductivity on 
the triangular lattice, there will still be a single Majo- 
rana bound state on the edge of the sample. An explicit 
example of the edge spectrum is shown in Fig. 8, where 
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(b): d+id-SC coexisting with 1x3 collinear SDW 




-3-2-10123 



k: momentum along 62 direction 

FIG. 9: (color online) The edge spectra with momentum k 
along 62 direction of a 2nd-NN d + id superconductor coex- 
isting with 1x3 collinear magnetic order with magnetization 
along z direction. Notations are the same as in Fig. 2 of 
the main text. Note the absence of Majorana bound state 
in the spectrum. The data are obtained for pairing energy 
A2 = 150 meV and magnetic mass m — 100 meV for the 
collinear order. 



0.01 




ft] -0.008 



k 2 : momentum along (1,1) direction 

FIG. 10: (color online) Edge spectra as a function of mo- 
mentum k 2 along the (1, 1) direction of the magnetic zone for 
a NN d x 2_ y 2 superconductor with (ir/2,ir/2) non-collinear 
magnetic order. The spectrum is zoomed in around k 2 — it 
where a flat band of zero energy states are localized on the 
two edges (blue and magenta), corresponding to the MBS. 
Red lines are the bulk states. Here we add a small 2nd NN 
imaginary pairing term Ak.i = 0.02i cos(fci). 



a single Majorana mode is localized on each of the two 
parallel edges and disperses across k = 0. 



C. Absence of Majorana fermions with collinear 
magnetic order 

As discussed in the main text for the generic case, a 
collinear magnetic order described by rl c i in Eq. (3) in 
the main text will not create a Majorana bound state in 
the vortex core or on the sample edge, although the bulk 
excitations are fully gapped. In Fig. 9, we show explic- 
itly the absence of Majorana mode in the edge spectrum 
obtained for a 2nd-NN d+ id superconductor coexisting 
with 1x3 collinear magnetic order with magnetization 
along z direction. 



D. Chirality of Majorana edge mode 

As discussed in the main text, when the pairing order 
parameter preserves time reversal symmetry the chirality 
of Majorana edge mode is C = sign[z • (v l x x Vy)] = 0, 
corresponding to a dispersionlcss MBS on the edge, as 
shown in FIG. 4 for the NN real d x 2_ y 2 pairing on square 
lattice. However in the presence of a magnetic order, 
time reversal symmetry is broken and a small imaginary 
pairing term Ai can be induced. This would result in a 
dispersing Majorana edge mode with C = ±1. Here we 
demonstrate this fact by calculating the edge spectrum of 
a d x 2_ y 2 pairing coexisting with (tt/2,it/2) non-collinear 
magnetic order, when a very small 2nd NN imaginary 



0.01 r 

+3 



o 




^ -o.ooe 



k 2 : momentum along (1,1) direction 

FIG. 11: (color online) Edge spectra as a function of mo- 
mentum k 2 along the (1, 1) direction of the magnetic zone for 
a NN d x 2_ y 2 superconductor with (tt/2,tt/2) non-collinear 
magnetic order. The spectrum is zoomed in around k 2 — n 
where a flat band of zero energy states are localized on the 
two edges (blue and magenta), corresponding to the MBS. 
Red lines are the bulk states. Here we add a small 2nd NN 
imaginary pairing term Ak,i = — 0.02i cos(fci). 



pairing A k ,i = ±0.02 i cos(fci) is induced. As shown in 
FIG. 10 and 11, the dispersionless MBS in FIG. 4 begins 
to disperse in the presence of the imaginary pairing term. 
When the sign of this imaginary pairing changes, the 
chirality of the Majorana edge mode changes accordingly. 
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